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Abstract 



Using a combination of renormalisation group and effective potential methods, we 
discuss the bound on the hghtest CP-even Higgs boson mass ruh in the next-to- 
minimal supersymmetric standard model. We find ruh < 146, 139, 149 GeV for nit = 
90, 140, 190 GeV. 
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1. Introduction 

1.1. Triviality 

The question of tlie origin of electroweak mass is one of the most urgent ques- 
tions of present day particle physics. The discovery of a particle which resembles 
the Higgs boson of the minimal standard model, and the measurement of its mass, 
will provide clues as to the nature of new physics beyond the standard model. In 
this report we shall be concerned with the question of how heavy the lightest neutral 
CP-even supersymmetric Higgs boson, h^, can be within the framework of super- 
symmetric grand unified theories (SUSY GUTs) In SUSY GUTs all the Yukawa 
couplings are constrained to remain perturbative in the region Msusy ~ 1 TeV to 
Mgut ~ 10^^ GeV. This constraint provides a maximum value at low energies for 
those Yukawa couplings which are not asymptotically-free, and is obtained from the 
renormalisation group (RG) equations together with the boundary conditions that 
the couplings become non-perturbative at Mqut ^ the so-called "triviality limit" . In 
the minimal supersymmetric standard model (MSSM) 0, the triviality limits provide 
a useful bound on the top quark mass mt. The upper bound on the mass, m/i, 
in the MSSM, including radiative corrections, has recently been the subject of much 
discussion^, ^, |^, |^, ||, p. However the MSSM is not the most general low energy 
manifestation of SUSY GUTs. It is possible that SUSY GUTs give rise to a low energy 
theory which contains an additional gauge singlet field, the so called next-to-minimal 



supersymmetric standard model (NMSSM) [0, |Tl|, 0. Here we shall concentrate 
on the question of the upper-bound on in the NMSSM which is obtained from 
triviality limits of Yukawa couplings. 

1.2. The NMSSM 

The NMSSM provides an elegant extension of MSSM which eliminates the 
problem, the appearance of a mass-scale in the superpotential. The NMSSM differs 
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from the MSSM by the presence of a gauge singlet field whose vacuum expectation 
value (vev) plays the role of the mass parameter fi in the MSSM, and is defined by 
the superpotential 



where the superfield = {tL^hL) contains the left-handed top and bottom quarks, 
and contains the charge conjugate of the right-handed top quark; Hi and H2 are the 
usual Higgs doublet superfields and N is the Higgs gauge singlet superfield; the ellipsis 
represents terms whose relatively small couplings play no role in our analysis (in 
particular, we assume that the bottom quark Yukawa coupling, h^, satisfies hi, <^ hf). 
The fields = (iJ°,/7f), i/J = {H 2,1^2) and N develop vevs which may be 
assumed to be of the form 



where z/i, 1/2 and x are real, and yz/f + z/f = z/ = 174 GeV. The low energy physical 
spectrum of the Higgs scalars consists of 3 CP-even neutral states, 2 CP-odd neutral 
states, and 2 charged scalars. A third CP-odd state is a Goldstone mode which 
becomes the longitudinal component of the , while a further two charged scalars 
become those of the W^^'s. 

An upper bound on the lightest neutral CP-even scalar in the NMSSM may 
be obtained from the real symmetric 3x3 neutral scalar mass squared matrix, by 
using the fact that the smallest eigenvalue of such a matrix must be smaller than the 
smallest eigenvalue of its upper 2x2 block. The resulting bound at tree-level is |r2| 



where tan/3 = ^, and A is regarded as a running parameter evaluated at Msusy- The 
upper bound on rrih is determined by the maximum value of \{Msusy), henceforth 




(1) 




< N >= X, 



(2) 



< Mz^ + {\^v^ - Mz^) sin^ 2(3. 



(3) 
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denoted Xmax- The value of Xmax is obtained by solving the SUSY RG equations for 
the Yukawa couplings ht, A and k in the region Msusy = 1 TeV to Mqut = 10^^ GeV 
P^ , P^ . If the Yukawa couplings ht, A and k are all initially large at Mqut then they 
approach low energy fixed point ratios |T^. However, if the boundary condition at 
Mgut is a fc, then larger values of X{Msusy) can be achieved. We have repeated 
the calculation of ref.]!^ and found that for ht^Msusv) = 0.5 — 1.0, Xmax = 0.87—0.70 



and for ht^Msusr) 1.06,Q Xmax (with A; = always). Radiative corrections 
to the tree-level bound in Eq.(2) have been considered in refs.[|l5, 16, 17|. In ref.|16 



these were estimated from a low energy RG analysis of the Higgs sector of the model 
between Msusy and a lower scale fi, assuming that only one Higgs boson has a mass 
below Msusy- 



In Section 2 |T3] we shall discuss a more general RG analysis in which both 



Higgs doublets and the Higgs singlet may be lighter than Msusy, making a simple 
approximation of hard decoupling below Msusy = ^TeV of the superpartners. In 
Section 3 ll9| we shall go on to consider the more general case in which the squarks 



have a more general mass spectrum, using the effective potential approach. Our 
numerical results for the bound are also presented in this section. Section 4 contains 
our concluding remarks. 



2. Renormalisation Group Approach ||18|| 

The basic assumption of this approach is that the squarks are degenerate at 
Msusy = 1 TeV while the Higgs bosons and top quarks have masses ~ 150 GeV. 
The effective theory below Msusy is just the standard model with two light Higgs 
doublets and a light Higgs singlet. Thus the Higgs potential at some low energy scale 

ht{MsuSY) < 1-06 is the triviality bound, which, together with rrit — ht{mt)vsm(3, where 
ht{mt) < 1.12, imphcs the bound nit < 195 GeV. 
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H < MsusY is given by the general expression 

- X4H2^H^\^ + X,\N\^\Hi\^ + \e\N\^\H2\^ 
+ X-r{N*^H^H2 + H.c.) + Xs\N\^ 

- mi{HiH2N + H.c.)-\m5{N^ + H.c.). (4) 

The running quartic couphngs Aj(/i) and the mass parameters mi{fi) must satisfy the 

following boundary conditions at Msusy 

1 1 
Ai = Aa = -(^2^ + ^1^), A3 = -(^2^ -^1^) 

A4 = A^ — '^92' ■, A5 = Ae = A^, A7 = — A/c, As = ■, 

mi = rriH^, m2 = rriHi, Tn3 = mN, 

7714 = XAx, m5 = kAk, (5) 

where Ax and Ak are soft parameters associated with the trilinear Higgs couplings in 
Eq.(l), and we have deferred a discussion of squark effects to Section 3. At energy 
scales /i below Msusy, the values of the quartic couplings may be obtained by solving 
the RG equations given in ref . |TB[ . The minimisation conditions implied by ^^g^'^'^" = 
and ^^q'^^" = allow us to eliminate the low energy parameters mi, m2, m^. The 
remaining masses and 777,5 are related to the parameters Ax and A^ at Msusy by 
Eq.(5). Below this scale we shall regard 7774 and 7775 as free parameters. 

The neutral CP-even (scalar) mass squared symmetric matrix in the basis 
1,2,3 = Hi,H2,N is 

I 2Xiul 2(A3 + A4)i^ii^2 2X^xi>i 

2(A3 + X4)ViV2 2X2^1 2XfiXV2 

\ 2X5x1^1 2X^X1^2 4A8X^ — 77750; 

/ tan/?[7774X — AyX^] — [7774X — AyX^] — ^ [7774a; — 2 AyX^^ 



+ 



— [7774X — X^x'^] cot /3[7774X — X^x'^] —^[m^x — 2A7X^] | (6) 
V - ^ [rriiX - 2 Atx^] - ^ [rmx - 2 Ayx^] ^ [m^x] 



The neutral CP-odd (pseudoscalar) mass squared symmetric matrix is 



[m4X — Xjx^] cot /9[m4X — Xjx"^] ^[Tn4X + 2\jx^] . (7) 

V flm^x + 2\7X^] ^[rriiX + 2\7x'^] ?,m^x\ ^[m^x - A^x"^] 



In ref . |]T8[ we obtained an upper bound on the mass of the hghtest neutral scalar 
Higgs boson h^, by using the fact that m\ must not exceed the lower eigenvalue of 
the upper 2x2 block matrix. The resulting upper bound is a complicated function 
of TTic^. It is easy to show that the bound reaches a maximum asymptotically for 



rric — > 00 [^. It is easy to show that the bound is maximised for A > Xmax and 
A; = 0, and so we shall use Xmax and A; = 0, as in the case of the tree level bound 
discussed previously. With this information in hand, a bound may be calculated from 
the values of 5\i = — \i{MsusY) obtained from numerically integrating the RG 
equations, decoupling the top quark below its mass and choosing /x = 150 GeV. For 
each value of m^, we have determined the value of ht{MsusY) and the corresponding 
value of Xmax which maximises the bound from a numerical analysis of the triviality 
condition as discussed previously. We defer a discussion of this bound until after 
squark effects have been considered. 

3. Squark Contributions [|T^ 

In order the calculate the shifts in the mass-squared matrices due to squarks, 
we employ the full one-loop effective potential. Much of our analysis of the squark 
spectrum is similar to that of the MSSM in Ref. P], whose notation we follow closely. 



A similar approach has also been followed by Ellwanger [^. The one-loop effective 
potential is given by 

ViiQ) = Vo{Q)+AVi{Q), 



where Vq{Q) is the tree-level potential at the arbitrary MS scale Q, Str denotes 



the usual supertrace, and Ai^ is the field-dependent mass-squared matrix. We shall 
take Q = Msusy = 1 TeV. B Vq is restricted to the pure Higgs part of the tree- 
level potential. Since the RG approach of Section 2 |TB| has already included the 
logarithmic effects of top quark and Higgs boson loops we exclude these particles 
from the supertrace. The logarithmic contributions to AVi{Q) from the top quark 
and Higgs bosons may be absorbed into Vo{Q) where Q = 1 TeV, to yield Vo(/u) where 
fi = 150 GeV, as in Eq.(4). It only remains to calculate AVi{Q), where Q = I TeV, 
involving squark contributions. 

To calculate the shifts in the Higgs mass matrices as a result of squark effects 
we require the field-dependent mass-squared matrices of the squarks. Taking the 
contribution to the soft SUSY breaking potential involving squarks to be 

AF.o/t = htAt{QH2t^ + h.c.) + ml\Q\^ + ml\i^\^ + ml\¥\^ , (9) 

where 6^ contains the charge conjugate of the right-handed bottom quark, and tildes 
denote the scalar components of the superfields, together with the superpotential in 
Eq.(l), leads to the field-dependent squark mass-squared matrix (ignoring contribu- 
tions proportional to gauge couplings and hi,) 

( ml + hj\H^\'^ XhtNH^ + htAtH^ -hfH^H^ \ 

2^ XhtNH'^ + htAHl ml + hl{\Hl\'' + HtH^) Xh^N Ht - htA^Ht 
-h^^H^H^ XhtNH{ - htAtH^ + h'^^H^H^ 

^ ml J 

(10) 

in the basis {t^, 6^, 6^} where a bar denotes complex conjugation. Manifestly, 
one eigenvalue is field-independent and may be discarded. Only the upper 2x2 
submatrix contributes to the CP-even and CP-odd mass-squared matrices, since the 
charged fields have zero vevs in order not to break QED, whereas the upper 3x3 
submatrix contributes to the charged mass-squared matrix. By differentiating AVi 

■^This value of Q is of course totally arbitrary. The value of 1 TeV is chosen to be consistent with 
our previous analysis. The results are independent of Q to one loop order. 
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once with respect to the fields we may obtain the shifts in the minimisation conditions 
induced by squark loops, and twice will yield the shifts in the mass matrices, i 

In the basis {Hi, H2, N} we have the following mass-squared matrices, after 
ensuring that the full one-loop potential is correctly minimised. The couplings Aj 
are those obtained from the potential in Eq.(4) renormalised at the scale n = 150 
GeV. For notational simplicity we drop the tildes. Moreover, we work in the basis of 
mass eigenstates {ti, ^2, ^2} (though, since hb, the bottom quark Yukawa coupling, 
is taken to be zero, 62 never contributes to the mass-squared matrices), where the 
mass of the ti squark is rrit^, and so on for the other squarks. The CP-even (scalar) 
mass-squared matrix is 

= + 6M^, (11) 

where is given in Eq.(6), and 



6M' 



/ A2 A2 A2 
' ^11 ^12 ^13 
.2 a2 a2 



\ A^ A^ 

\ ^13 ^23 ^ 



33 




A?2 Al, Al, + -1 cot/3 A^. (12) 



The Ajj and A^ are given by 
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Al = ^^hl{Xx).A.f{mlmlj, 



All = ^.htuliXxY. ( ^I.ti'^^g^ ) 9irnl<l 



A^ 
^22 



3 4 2 A mlml ^ 2At{At + Xx cot p) ^ ml' 
- Svr^^*"- rM|^,,+ ml -ml 

+ ^/^.v,-.f ^-^l-2^^1f^^ V^K.<). 



— m; 



2 



a 2 3 4 2 /\ n2 ( At + Xx cot (3 \ 2 2\ 



Stt^ Y ml — ml J \ ml ml — m 



■ti 



§In fact, this is only approximately true due to Higgs self-energy corrections; these are expected 
to be small for the lightest Higgs bosons fl. 
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^13 



^23 



m| - ml 



3 ,4 2/\ wx . (At + XX cot P\ 2 2> 

3 
8^ 

'^l A2 



h].{\x).{\vi).f{m\,m\] 



X 



^125 



where the functions / and 5^ are defined by 



t2J 



9(^1, ml) 



ml - ml 



ml - ml 



2 1 ^^ti 2 2 1 

my- m — m, — m, m 



m 



/\/f2 



i^ + m^ 



t2 



m 



(13) 



(14) 



If the squarks are degenerate in mass, then Eq.(lO) imphes that At + Ax cot /5 = 
0. Furthermore, if mt-^ = mt^ = Msusy, then Eqs.( 12-14) imply that the squark 
contribution to the CP-even mass-squared matrix vanishes. This is the hmit of our 
previous analysis in Section 2 |T8[]. 



The CP-odd (pseudoscalar) mass-squared matrix is 



where is given in Eq.(7) and 



6M^ 



tan/? 1 f 
1 cot/5 ^ 



1^2 



1^1 



,.2 



A^. 

p 



Finally, the charged mass-squared matrix is 
/ tan j3 1 



1 



cot /3 



{m^x — Xjx — A4Z/ii/2 + A^), 



where 



and 



167r2 



E 



m„ In 



m„ 



d'^m^ 



ma<^{mti,mt2,mt-^} 



/\/f2 



- 1 



OH^dH; 



-)- I vevs 1 



d'^ml 
dH^dHt 

d'^ml 
dH^dHt 

d'^ml 
dH{dHt 



hjpl^Xxf.cotp hl{Xx).At 



2 ' 



{ml - ml) {ml - mfj ml - ml 
hjvl{\xf.coiP hl{\x).At 



+ 



{ml - ml ) {ml - ml ) ml - ml ' 
hlul.{\xf.cot[3 



{ml -ml){ml - ml)' 



(15) 



(16) 



(17) 



(19) 
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The bound on the hghtest CP-even Higgs mass is a consequence of the fact 
that the minimum eigenvalue of M"^ is bounded by the minimum eigenvalue of the 
upper 2x2 submatrix of M^. Using Eq. (|T7p to obtain the physical charged Higgs 
mass squared, m^, in terms of 7/14, 

rriiX - = -{ml + X^u^) sin 2/3 - A^, (20) 

we may eliminate mi from in favour of m^, and thus write the upper 2x2 
submatrix in the form 

= M'2 + 6M'\ (21) 

where 



M'2- ( 2(A3 + A4)z/iz/2 \ / tan^ Vr^^ , , ,,2n „• 

(22) 



and 



It is clear from the form of the second term on the right-hand-side of Eq. (|23|) that 
the factor — A^ does not change the bound at all, but merely serves, for a fixed 
tan/3, to shift the bound, when plotted as a function of mc, to the left or right. We 
may then drop this term, since its presence may be re-parametrised by a shift in the 
free parameter mc- In this way m^^ is eliminated. 

From Eq. (^^ and Eq. (P3| ) it is a simple matter to determine the shift in the 
minimum eigenvalue of M'^ due to squark effects. Defining A = (M'^)ii, B = (M'^)i2 
and C = (M'^)22, the shift is given by 



.2 ^ .2 (A-C)(A?,-Ay +4i?A?2 
^11 + ^22 



(24) 



'{A-Cy + AB^ 

Because this shift is a one loop effect, the couplings Aj in A, B and C may be evaluated 
at any renormalisation point, since the difference between a coupling evaluated at two 



different renormalisation points is also a one loop effect, thus giving an overall error 
at the two loop level, which we neglect in our approximation. 



In our previous analysis we used triviality limits on the couplings ht, X 



and k [^3[ |T4| to determine the bound on the lightest CP-even Higgs boson mass 
m° not including general squark effects, that is, with 5M'^ = 0. Let the triviality 
limit on A for a given hf be Xmax- In table |l]we reproduce the values of ht and Xmax 
which generated the bound m° of our previous analysis — all other values, for a given 
top quark mass, resulted in smaller values of the bound. It transpires that = in 
all cases. S In order to maximise the bound resulting from Eq. p^ ) we take the 
limit Trie = oo, thus eliminating TfJiQ clS db parameter. As a function of rric the bound 
approaches its maximum asymptotically as oo. The approach is rapid, with 

rric ~ 200 GeV being a good approximation to = oo. Thus, the bound resulting 



from the purely formal procedure of taking the limit rric = oo is not unphysical [ITH 



Had this not been the case, then certainly we would have an upper bound, but one 
perhaps not capable of realisation in an actual spectrum. 

We use the values of ht and Xmax in table |l| to calculate the shift in the 
bound on the lightest CP-even Higgs boson mass-squared given by Eq. (p^). The 
shift in the bound on the mass is denoted by 6mh. This is a hideously complicated 
function of many parameters, including the squark masses m^^ and m^j, the soft 
SUSY breaking parameter At, the vev of the gauge singlet field x, and rric- As a 
function of rric the shift appears, numerically, to be maximised for rric small, typically 
between 100 GeV and 200 GeV, with the difference between the shifts at these two 
points being less than 1 GeV. Thus, we set rric = 200 GeV in the general squark 
contributions to the bound. Doing this enables us to retain the bound from our 
previous analysis derived from Eq. (^21) in the limit rric = oo, and simply maximise the 



^Taking A: = gives rise to an axion in the physical spectrum. However, a small, non-zero value 
of k is sufficient to give the would-be axion a mass of several tens of GeV. Such a value of k does 
not invalidate our calculations. 
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Table 1: Lightest CP-even Higgs mass bound in the NMSSM including general squark effects. In 
row 1 is the top mass, rrif, in GeV; in rows 2 and 3 those values of h^i^MgusY) and XmaxiMgusy): 
which produce the bound in our previous analysis, this being in row 4, in GeV; in row 5 is the 
contribution to the bound resulting from a general squark mass spectrum rather than degenerate 
squarks with mass MsuSYj in GeV; in row 6 is the lightest CP-even Higgs mass bound including 
general squark effects, in GeV; in rows 7 and 8 are those value of nit^ and At, in GeV, which 
generate the bound in row 6; = 1 TeV. 



mt 


90 


100 


110 


120 


130 


140 


150 


160 


170 


180 


190 


ht 


0.60 


0.67 


0.73 


0.79 


0.84 


0.88 


0.92 


0.95 


0.98 


1.00 


1.03 


^max 


0.87 


0.85 


0.83 


0.81 


0.79 


0.77 


0.74 


0.71 


0.67 


0.63 


0.50 


ml 


145 


143 


140 


137 


135 


132 


129 


126 


124 


123 


126 


Sruh 


1 


2 


3 


4 


6 


8 


10 


13 


16 


20 


23 


nih 


146 


145 


143 


141 


140 


139 


139 


139 


140 


143 


149 


mti 


800 


770 


750 


720 


700 


670 


650 


630 


610 


590 


570 


At 


940 


1000 


1040 


1130 


1160 


1270 


1330 


1410 


1500 


1630 


1820 



squark contributions separately and add them to our previous bound to yield the new 
bound rrih = + Srrih- The x dependence is not too strong, with the ratio r = x/u 
typically taking values close to 10 in order to maximise the bound. Where this is not 
the case, the difference between the bound at its maximum, as a function of x, and 
that at r = 10 is less than 1 GeV. Thus, we set r = 10, this giving a reliable indication 
of the maximum shift. We calculate the shift for various ranges of squark masses and 
a range of values of At, and record the maximum value obtained. The squark masses 
are allowed to vary between 20 GeV and 1000 GeV in steps of 10 GeV. (1000 GeV is 
the upper limit since we take Msusv = 1000 GeV.) Given this restriction, the value 
of At which maximises the general squark contribution to the bound never exceeds 2 
TeV. We impose the constraint 2htV2{At + Ax cot (3) < |m^^ — m^^ | which follows from 
the form of the squark mass matrix in Eq. ([T0|). 

In table |l] we show the bound on the lightest CP-even Higgs mass mh including 
general squark effects (and m° for comparison). Table |l] also shows the values of m^^ 
and At used to generate the bound. As mt increases, mt^ monotonically decreases, 
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while At monotonically increases. The other squark mass, rrit^, takes the value 1 
TeV over the whole range of top quark masses — that one squark mass takes on 
its maximum permitted value to maximise the general squark contribution to the 
bound can be seen analytically. It will be noticed that these squark masses and 
At are somewhat larger than typically expected from GUT scenarios. However, we 
feel it wise not to restrict ourselves too excessively to particular prejudices regarding 
physics as yet unknown (though, of course, the notion of triviality does require the 
assumption of a SUSY desert up to the unification scale). 

Our calculations suggest that the universal upper bound on the lightest CP- 
even Higgs mass is 149 GeV. Table reveals that for large values of rrit squark effects 
may contribute up to 6mh = 23 GeV, but for small rrit the squark effects are small, 
as expected. We emphasise that we have not re-maximised the complete one loop 
corrected bound, but rather made the assumption that the values of ht and Xmax 
which generated the bound of our previous analysis are not significantly modified by 
the inclusion of a general spectrum of squark masses. 



4. Conclusions 



We have combined an RG analysis of the Higgs sector of the NMSSM |1 



with a calculation of general squark effects |19], and determined a bound on the mass 



of the lightest CP-even state of mn < 146, 139, 149 GeV for = 90, 140, 190 GeV. 
Our calculations indicate that the effects of a general squark spectrum can be very 
significant for large top quark masses. One reason why, for large top quark masses, 
squark effects are important is that there exist finite, one-loop diagrams with vertices 
containing the uncontrolled soft SUSY breaking parameter At. These diagrams can 
give large contributions to the Higgs boson masses. We close by indicating that 
in the NMSSM there are similar diagrams involving loops of Higgs bosons, again 
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with vertices containing soft SUSY breaking parameters, this time Ax and Ak- The 
existence of these diagrams is directly due to the gauge singlet field N; they do not 
occur in the MSSM. We see no reason in principle why these new diagrams should 
not give rise to similarly large effects, except for obvious factors of 3 due to colour. 
We are currently in the process of estimating these effects [^. 
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